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ABSTRACT 

ECG is a non-invasive technique is used as a primary diagnostic tool for cardiovascular disease. The 
proposed method will help to acquire an ECG signal by generating a synthetic ECG from a patient and to 
analyze and detect the anomalies which are also called as abnormalities. Abstract Data obtained from 
electrocardiogram (ECG) signals provides invaluable tools for diagnosing cardiac disorders. However, ECG 
signals recorded from electrocardiograph are usually corrupted by noise attributed to several factors. Finite 
Impulse Response Filtering Techniques will eliminate noise in the signal to some extent,but te resultant output 
signal has a noise component. To help solve these problems, we simulated a simple but inexpensive and easy-to- 
implement MATLAB model that generates ECG signals and gives us mathematical control over the ECG signal. 
Our model fuses Mathematical functions in MATLAB with physiological data. 
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I. INTRODUCTION 

Electrocardiogram (ECG) is a diagnosis tool that reported the electrical activity of heart recorded by 
skin electrode. The morphology and heart rate reflects the cardiac health of human heart beat. It is a non- 
invasive technique that means this signal is measured on the surface of human body, which is used in 
identification of the heart diseases. Any disorder of heart rate or rhythm, or change in the morphological pattern, 
is an indication of cardiac arrhythmia, which could be detected by analysis of the recorded ECG waveform. The 
amplitude and duration of the P-QRS-T wave contains useful information about the nature of disease afflicting 
the heart. The electrical wave is due to depolarization and re polarization of Na + and k ions in the blood .The 
ECG signal provides the following information of a human heart. 

A typical ECG signal shows the oscillations between cardiac contractions (systole) and relaxations (diastole) 
states as reflected in a heart rate (HR). Thus the ECG signal determines the number of heart beats per minute. 

A number of important events characterize cardiac functions. Atrial and ventricular depolarization/re - 
polarization takes place for each heart beat. The cardiac cycle is associated with portions of the heart becoming 
positively charged, while the remaining parts become negatively charged interchangeably. This potential 
difference generated initiates the flow of current. A typical ECG signal depicts a series of waveforms which 
occur in a repetitive order. The waveforms are initiated from the isometric line, from which a deflection 
indicates electrical activity. One normal heart beat is represented by a set of three recognizable waveforms that 
start with the P-wave, followed by the QRS complex and ends with the T-wave. The relatively small P-wave is 
initiated by the depolarization of the atrial muscles and is related to their contraction. 

The large QRS -wave complex, made up of three waves, is caused by the depolarization of the 
ventricles and is connected to their contraction. Atrial re -polarization happens during the depolarization of the 
ventricles but its weak signal is undetected on an ECG. The T-wave is caused by currents flowing during the 
repolarization of the ventricles. A normal cardiac cycle of an individual at rest consisting of all waveforms 
(form P -T waves) spans 0.8 seconds. 

ECG signals are generated using electrocardiographs. Such signals are usually vitiated by several 
sources of noise which include (i) electrical interference from surrounding equipment (e.g. effect of the 
electrical mains supply), (ii) measurement (or electrode contact) noise, (iii) electromyographic (muscle 
contraction), (iv) movement artifacts, (v) baseline drift and respiratory artifacts and (vi) instrumentation noise 
(such as artifacts from the analogue to digital conversion process). One method of dealing with corrupt ECG 
signals is through the use of signal filtration systems. 
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II. METHODOLOGIES AND MATHEMATICAL EQUATIONS 

In order to eliminate the perennial problem of dealing with noise from the tainted ECG signals, 
dynamic functions are usually developed for the generation of artificial synthetic ECG signals. Such functions 
are usually equation based. The work of McSharry and Clifford relies on three coupled differential equations to 
generate ECG signals. We summarize their method below. McSharry and Clifford captured the spectral 
characteristics of beat-to-beat RR intervals or RR tachogram, including both the oscillation in the RR tachogram 
resulting from parasympathetic activity, in synchrony with respiration (Respiratory Sinus Arrhythmia), and the 
waves in arterial blood pressure (Mayer waves), using a bi-modal spectrum made up of the sum of two Gaussian 
functions. 
This is given by: 

Where fl and f2 are the two means, and cl and c2 are the corresponding standard deviations. In spectral analysis 
of the RR tachogram, two critical frequency bands, usually referred to as the low -frequency (LF) band (0.04 to 
0.15 Hz) and high-frequency (HF) band (0.15 to 0.4 Hz) are considered. 

The power in the LF and HF bands are denoted by o{ and erf respectively and the variance is represented by 

rJ 2 = rJ^ + of 2 



Consequently, the LF/HF ratio is given by ffj/ffj 3 

The function S (f) gives the spectrum of a time series denoted by T (t). In order to obtain T (t), the 
inverse Fourier transform is applied on a sequence of complex numbers whose amplitudes are given by *jS(f) 
, such that the phases are randomly distributed between the interval 0 and 2n. Next, an appropriate scaling 
constant is chosen to multiply the resulting time series and then an offset value is added. This makes it possible 
to assign any required mean and variance to the resulting time series T (t) based on the initial series. 

It also specifies different realizations of the random phases by simply changing the seed of the random 
number generator. The resulting series will inherit the same temporal and spectral properties of the original. To 
complete the process, a dynamical model is then used to obtain the quasi-periodicity property of the ECG signal. 
This is guaranteed by making sure that the model has an attracting limit cycle. One complete revolution around 
the limit cycle in the x-y plane mimics a heart-beat. The PQRST peaks and troughs of the ECG signal are 
captured by a series of exponential functions inherent in the model, while the PQRST are the extremes of the 
signal are specified by five angles correspondingly denoted by 0p, 0q, 0r, 0s, and 0t. The mathematical model is 
represented by the following system of differential equations. 

x = ax-coy 4 

y = ay + cox 5 

z = - EiEtP,Qjw;rj Aeiex P (~ Ae2 1 /2b2 1) ~ ( z_z * ) 5 

Where, a =1-V(x2 -y2), A0i = (0 -0i) mod (2n) 9 0 = atan2(y, x) is the four quadrant arctangents of the elements 
of x and y, ranging over [-71, 7i] and co is the angular frequency of the trajectory in its motion around the limit 
cycle and is related to the heart beat rate as 27if [9]. MATLAB Is used to develop and test "SynECG", the 
artificial ECG generating function. 

The five inbuilt MATLAB functions used to develop "SyECG" are "sgolayfilt", "kron" "ones", 
"round", and "linspace". These functions are described briefly below. Function sgolayfilt () is known as the 
Savitzky-Golay Filter, y = sgolayfilt(x, k, f) applies a Savitzky-Golay FIR smoothing filter to the data in vector 
x. If x is a matrix, sgolayfilt operates on each column. The polynomial order k must be less than the frame size, 
f, which must be odd. If k = f -1, the filter produces no smoothing. Savitzky-Golay smoothing filters (also called 
digital smoothing polynomial filters or least-squares smoothing filters) are typically used to "smooth out" a 
noisy signal whose frequency span (without noise) is large. It also helps preserve the peaks and valleys of ECG 
signals better than a standard FIR filter. 

III. SIMULATION RESULTS 

We use the SynECG to generate 10-seconds synthetic ECGs for a number of heart rates of a resting 
individual with a peak voltage of 1.2mV. Figure 1.2 Illustrates how the function is summoned in the MATLAB 
command window; Figures 2 through 4 show the artificial ECGs generated for 40bpm, 50bpm, 60bpm, 70bpm, 
80bpm, 90bpm, lOObpm and 1 lObpm heart rates at a desired peak voltage of 1.2mV. 
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We note that heart beats in a regular rhythm is usually between 60 and 100, that is when the signal-to- 
noise ratio is usually quite good in a person at rest. Basically, the ECG is a piece wise continuous graph of 
potential difference (in mV) against time (in seconds). The Synthetic ECG generated for 95bpm displays 
approximately 17 cardiac cycles in 10 seconds. 

With the SynECG, we have been able to generate a 10-seconds artificial ECG for some selected heart 
rates. The 60pbm generated 10-seconds ECG displays an approximate 10 cardiac cycles. 

As heart rate increases, the number of approximate cardiac cycles in the 10-seconds ECG also 
increases. Variation in the heart rates time intervals can be monitored with this increase in heart rate. A careful 
observation of the 60bpm ECG generated by SynECG reveals that one cardiac cycle takes 0.8 seconds as human 
physiology suggests. We observe from the generated ECGs that as heart rate increases; the time taken for a 
cardiac cycle reduces. Similar observations were made for 80bpm and 95bpm. The various waveforms which 
constitute a cardiac cycle vary in length as heart rates change from time to time for the resting individual. 

Artificial Electrocardiogram 
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Figure 2: Synthetic ECG signal generated from SYNECG using 40bpm. This displays approximately 10 
cardiac cycles in ten seconds. 
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Artificial Electrocardiogram 




Figure 2: Synthetic ECG signal generated from SYNECG using 50bpm. This displays approximately 10 
cardiac cycles in ten seconds. 
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Figure 3: Synthetic ECG signal generated from SYNECG for 60bpm. This displays approximately 13 
cardiac cycles in ten seconds. 
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Figure 4: Synthetic ECG signal generated from SYNECG using 70bpm. The graph displays 
approximately 17 cardiac cycles in ten seconds. 
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Figure 5: Synthetic ECG signal generated from SYNECG for 80bpm. This displays approximately 13 
cardiac cycles in ten seconds. 
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Figure 6: Synthetic ECG signal generated from SYNECG for 90bpm. This displays approximately 13 
cardiac cycles in ten seconds. 
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Figure 7: Synthetic ECG signal generated from SYNECG for lOObpm. This displays approximately 13 
cardiac cycles in ten seconds. 
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IV. CONCLUSION 

Finally we have also worked with the SynECG, with which we have been able to generate a 10-seconds 
artificial ECG for some selected heart rates. The 50bpm generated 10-seconds ECG displays an approximate 10 
cardiac cycles. As heart rate increases, the number of approximate cardiac cycles in the 10-seconds ECG also 
increases. Variation in the heart rates time intervals can be monitored with this increase in heart rate. A careful 
observation of the 60bpm ECG generated by SynECG reveals that one cardiac cycle takes 0.8 seconds as human 
physiology suggests. We observe from the generated ECGs that as heart rate increases; the time taken for a 
cardiac cycle reduces. Similar observations were made for 80bpm and 95bpm. The various waveforms which 
constitute a cardiac cycle vary in length as heart rates change from time to time for the resting individual. 

With the SynECG function, synthetic ECG signals of various heart rates can be generated at a desired 
peak voltage for off-line study of the electrical activity of the heart. Our function made use of the function 
"sgolay" which is specifically geared towards ECG noise reduction and uses Savitzy-Golay polynomial filtering 
which is well-suited to smoothing of ECG data and preserve the peaks and valleys of the ECG signals better 
than a standard FIR filter. We have been able to show that our function generates synthetic ECG for 1 0 seconds 
by inputting heart rate (per minute) and the desired peak in millivolts. 
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